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ABSTRACT 



We show that the evolutionary track of a low-mass red giant should make an 
P^ ■ extended zigzag on the Hertzsprung-Russel diagram just after the bump lumi- 

nosity, if fast internal rotation and enhanced extra mixing in the radiative zone 
bring the temperature gradient close to the adiabatic one. This can explain both 
the location and peculiar surface chemi cal composition of Li-rich K giants studied 



O. , , 

^ . by iKumar. Reddy. fc Lambert I ( 120111 ) . We also discuss a striking resemblance 

j^ . between the photometric and composition peculiarities of these stars and giant 

components of RS CVn binaries. We demonstrate that the observationally con- 
strained values of the temperature gradient in the Li-rich K giants agree with 
QQ . the required rate of extra mixing only if the turbulence which is believed to be 

^ . responsible for this extra mixing is highly anisotropic, with its associated trans- 

port coefficients in the horizontal direction strongly dominating over those in the 

cn . vertical direction. 
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C^ . 1. Introduction 

During their first ascent along the red giant branch (RGB), low-mass stars are known to 
experience extra mixing in their convectively stable radiative zones, between the hydrogen 
burning shell (HBS) and the bottorn of t he convective envelope (BCE). This occurs above 



the bump luminosity ( jGratton et al.ll2000l ). after the HBS has crossed and erased a chemical 
composition discontinuity left behind by the BCE at the end of the first dredge- up (FDU). 
As a result, the atmospheric ^^C/^^C ratio and, in metal-poor stars, carbon abundance 
resume their post-FDU declines with increasing luminosity. Extra mixing is needed to explain 
the observed low carbon isot opic ratios in evolved stars with initial masses up to 2.2 Mq 



( ICharbonnel &: Lagardell2010l ). Another abundance anomaly that has also been associated 
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with the start of RGB extra mixing at the bump luminosity i s the hthium enrichment that 
has b een detected in a small fraction (1-2%) of K giants (ICharbonriel &: Balachandran 



2000l ). To produce ^Li via the ^Be mechanism (jCameron fc Fowlerlll97ll ). extra mixing has 
to be nearly three orders of n aagnitude faster in these stars than in other low-mass red gian ts 
above the bump luminosity ( ISackmann. fc Boothroyd Ill999l : iDenissenkov. fc Weiss 1120001 ) . 



The problem of the origin of Li enrich ment in K giants has recently been addressed 
again by iKumar. Reddy. &: Lambert I (1201 ll . hereafter KRL). They observed 2000 K giants 
with determined Hipparcos parallaxes and found that 15 of them are Li-rich. However, the 
most interesting new result is that many of these Li-rich K giants are located well below 
the bump luminosity, close to the red horizontal-branch clump region of the Hertzsprung- 
Russell diagram (HRD), where the low- mass stars should arrive much later, after they will 
have experienced the He-core flash. The newly discovered Li-rich K giants also exhibit very 
low ^^C/^'^C ratios, approaching the equilibrium value for the CN cycle in four of them. 
Analysing these data, KRL have proposed that the lithium enrichment in these K giants 
occurred during the He-core flash rather than at the bump luminosity. 

In this Letter, we present an alternative explanation; namely, we assume that these 
stars have already reached the bump luminosity and, therefore, that RGB extra mixing is 
now operating in their radiative zones. However, for some reason, presumably because of 
their very rapid internal rotation, this mixing is so efficient that not only can it trigger 
the ^Be mechanism, but its driving force or associated heat transport can also modify the 
radiative zone's thermal structure, resulting in a temperature gradient that is closer to the 
adiabatic one. In this case, a bump-luminosity star in which such enhanced extra mixing 
has just started should make an extended zigzag on the HRD towards much lower lumi- 
nosities, comparable to those of red clump stars, before resuming its asce nt along the RGB 
( IDenissenkov. Chaboyer. fc Li II2006I : IDenissenkov. fc Pinsonneaultll2008bl ). This hypothesis 
is supported by a striking resemblance between the photometric and chemical composition 
peculiarities of the primary (red-giant) components of the RS CVn binaries and the Li-rich 
K-giants studied by KRL. 



2. Photometric Peculiarities Caused by Enhanced Extra Mixing 

The RS CVn binaries are close stellar systems in which the primary components are 
low-mass red giants. In many cases, the RS CVn giant components have synchronized their 
spin and orbital periods and, as a result of their fast rotation caused by the tidal sp i nning 



up, they show signa t ures of chromosphe ric act ivitv (IFekel et al. 



Fekel fc Henrvl l2005l ). Morel et al.l (120041 ) . and IDenissenkov. Chabover. fc Li I (l2006l ) have 
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noted a remarkable photometric peculiarity intrinsic to the RS CVn giants and other K- 
giants in tidally locked binaries — most of them reside on the lower RGB, below the bump 
luminosity. 



Denissenkov. fc Pinsonneaultl ( l2008bl ) have shown that, after having reached the bump 



luminosity, a low-mass star will make a long excursion towards much lower luminosities if 
the temperature gradient V = {d\nT/d\nP) in its radiative zone increases and takes a 
value between the radiative and adiabatic temperature gradients, Vrad < V < Vad- Such 
a deviation of V from Vrad can be explained in two ways. On the one hand, a rotational 
distortion of level surfaces in the radiative zone leads to an increase of V = Vrad to V ~ 
(1 + 2e)Vrad, where e ~ {fl'^r^/SGMr) < ^crit is the ratio of the equatorial centrifugal 
acceleration to gravity divided by 3, and £crit = 0. 24 is its critical value at which the R oche 
equipotential surface breaks up at the equator (iDenissenkov. &: VandenBerg I l2003bl ) . It 
is this effect that was proposed by IDenissenkov. Chaboyer. fc Li I (120061 ) to be responsible 
for the photometric peculiarity of red giants in tidally locked binaries. On the other hand, 
fast differential rotation and its related hydrodynamic or magnetohydrodynamic instabilities 



can be driv ing mechanisms for extra mixing in stellar radiative zones (e.g., IZahru Il992 



Spruitlll999l ). Besides chemical elements, this extra mixing should also transport heat. The 
efficiency of the latter process in the standard mixing length theory is estimated by the 
quantity F ^ 7(Dmix/-ft'), where -Drnix is the rate (diffusion coefficient) of extra mixing, K is 
the thermal diffusivity, and 7 is the factor determined by the geometry of the fluid elements. 
In this case, the temperature gradient lies between Vrad and Vad, and its value is given by 

V ^ (1 - /)Vrad + /Vad, whcre / = 6rV(i + r + er^) dMaedeJIiggsl ). 



The giant components of the RS CVn binaries selected by us from the paper of iMorel et al 



( I2OO4J ) for comparison and the Li-rich K giants discovered by KRL both have nearly solar 
metallicities and masses close to 2Mq. Therefore, we compare their positions on the HRD 
(the red and blue star symbols in Fig. [T]) with the evolution of two stellar models calcu- 
lateqj for the solar heavy-element and helium mass fractions, Z = 0.019 and Y = 0.28, 
and initial masses 1.7 Mq and 2.2 Mq (the black and green curves in the same flgure). The 
neighbouring open symbols along the evolutionary tracks are separated by one million years. 
Their concentrations on the green curve at low T^s are increased at the bump luminosity, 
logiQ{L/LQ) ~ 2.25, and at the red clump, logiQ^L/LQ) ^ 1.7. Before the star has evolved 
beyond the bump luminosity, its radiative zone contains a strong gradient of the mean molec- 
ular weight, V^ = dlnfi/dlnP > 0, that is believed to prevent any extra mixing. The blue 
star symbols are all located below the bump luminosity, in a region of the HRD that the 



-'^AU the stellar evolution computations for this paper have been done with the MESA code (|Paxton et al. 



20111 ) that is freely available at http://mesa.sourceforge.net 
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standard evolutionary track (the green curve) of a low-mass red giant with a chemically 
uniform radiative zone can reach only after the He-core flash. This explains why KRL have 
inferred that both the production of Li and the reduction of the carbon isotopic ratio in their 
Li-rich K-giants occurred during the He-core flash. 

However, the ^Be mechanism of surface Li enrichment requires strongly enhanced extra 
mixing with a rate (diffusion coefficient) increased by nearly three orders of magnitude 
compared to that of c anonical extra mixing operating in t he majority of RGB stars above 



the bump luminosity ( jDenissenkov. fc VandenBerg Il2003al ). We surmise that such mixing 



or its associated fast rotation should modify V in the radiative zone bringing its value closer 
to Vad- The red curve originating from logiQ{L/LQ) ^ 2.25 in Fig. [1] shows a bifurcation of 
the evolution of the 2.2 Mq star after the bump luminosity caused by our assumption that 
the start of enhanced extra mixing with -Dmix = 2 x 10^^ cm^ s^^ in this star has initiated a 
convective heat transport in its radiative zone with the efficiency F = 0.75. It is seen that 
the star spends quite a long time at luminosities close to those of the Li-rich K-giants. We 
have assumed that the enhanced extra mixing has the same maximum depth, rmix = 0.05 Rq, 



as the canonical one (JDenissenkov. fc Pinsonneaultll2008al ). Before the HBS, advancing in 



mass outwards, has crossed the chemical composition discontinuity left behind by the BCE 
at the end of the FDU, the mixing is allowed to operate only in the chemically uniform 
zone between the BCE and the current location of the discontinuity at r > rmix- Note that 
it is the increase of V, not the effect of mixing, that forces the star to make the extended 
zigzag. Fig. [T] also shows that the region occupied by the majority of tidally locked giants 
in close binaries, including the RS CVn giant components (red star symbols), is reached by 
the evolutionary U-turn of the 1.7 Mq star (red curve) when F = L 



3. Chemical Composition Peculiarities Caused by Enhanced Extra Mixing 

The parameter values -Dmix = 2 x lO^^cm^s^^, rmix = 0.05 i?©, and F = 0.75 lead 
to both Li enrichment^ and reduction of the ^^C/^'^C ratio consistent with those reported 
by KRL for the Li- rich K giants (solid red curves in Fig. [2] and Fig. [3]). Furthermore, an 
increase of the efficiency of convective heat transport to F = 1 shifts the HRD location of 
the Li enrichment and, especially, that of the ^^C/^^C sharp decline to the edges of their 
corresponding observational domains (dashed red curves in the same figures). Therefore, if 
our hypothesis is correct, F should be close to 0.75 and not exceed 1 by much. 



^For the Li abundance, we use the notation £(^Li) = log]^Q[n(^Li)/n(H)] + 12, where n is the number 
density of the corresponding nucleus. 
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A similar value of F is needed to explain the photometric peculi arity of red giants in 

tidally locked binaries, including the RS CVn giant components (Fig.[T]). iDenissenkov. Chaboyer. &: Li 
( 120061 ) have predicted that some of these stars, those that have already reached the bump lu- 
minosity and are now making extended zigzags towards lower luminosities as a result of their 
fast rotation, should have carbon isotopic ratios smaller than the st andard post-FPU value of 
i2q^i3q ;-^ 25. Recently , two such stars have actually been found (JTautvaisiene et al.ll2010l : 



Barisevicius et al.ll2010l ): the chromospherically active RS CVn-type stars A And (HD 222107) 
and 29 Draconis (HD 160538) are located below the bump luminosity but have ^^C/^'^C = 14 
and 16, respectively. 

When comparing the Li-rich K giants with the RS CVn binaries, it is important to bear 
in mind that a hunt after the former is biased towards finding the stars that have already 
reached the bump luminosity and are now experiencing enhanced extra mixing, while a 
search for the latter selects both the pre- and post-bump luminosity stars. Therefore, it is not 
surprising that the third RS CV n-type star, 33 Piscium, studied by the same group does not 
show signatures of extra mixing (JBarisevicius et al.ll201ll ). It would be interesting to measure 
carbon isotopic ratios in the following RS CVn stars: HD 19754, HD 182776, HD 202134, 
HD 204128, and HD 205249. They have relatively high (for red giants) projected rotational 
velocities (vsi ni > 7 km s~^) and probably strongly increased atmospheric abundances of 



(IDenissenkov. fc VandenBerg 



Na and Al (Morel et al.l 12004 ) , whic h might be a manifestation of enhanced extra mixing 



2003ah . 



For comparison, green curves in Fig. [2] and Fig. [3] show the evolutionary changes of the Li 
abundance and ^^C/^'^C ratio at the surface of the 2.2 Mq star produced by the FDU (before 
the bump luminosity) and thermohaline convectio n driven by '^He burning t hat reduces fi 
in the tail of the HBS after the bump luminosity (ICharbonnel fc Zahn 20071). Th e rate of 



the latter is approximated by D^ = 27r^a^[V^/(Vrad — Vad)]-^' (lDenissenkovll2010al ) with the 



observationally constrained value of the aspect ratio, a = l/d = 10, of a fluid element, where 
/ and d are its length and diameter. Note that direct numerical simulations of the ^He-driven 
thermohaline convection in a bump-luminosity red giant give an estimate of a ~ 1, which 
does not support the hypothesis that it is the princ i pal mechanism of RGB extra mixiri g 
(JDenissenkovl l2010al : IDenissenkov. &: Merryfieldl l201ll : iTraxler. Garaud. &: Stellmach I l201ll ) . 



4. Discussion 



For the majority of low-mass red giants, that are experiencing extra mixing above the 
bump luminosity, the observed evolutionary changes of their surface chemical composition 
are reproduced reasonably well either with the above-mentioned thermohaline diffusion coef- 
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ficient D^ that uses a = 10, or with the diffusion coe fficient Dm\^ ^ O.Ol/f-O.liT and app ro- 
priately chosen mixing depth, e.g. rmjx = 0.05 i?© ( iDenissenkov. &: Pinsonneaultll2008al Jbl). 
For the radiative zone of our 2.2 Mq bump-luminosity model, these diffusion coefficients are 
plotted in Fig. H] (the green and black curves). Their corresponding efficiencies of heat trans- 
port are negligibly small. Therefore, in the majority of cases, the (canonical) RGB extra 
mixing should not lead to noticeable changes in the evolution of these stars on the HRD. 

On the contrary, the value of -Dmix = 2 x 10^^ cm^s~^, that is required to understand 
the origin of Li- rich K giants in our models, results in F = ■j^D^i^/K) :» 1, at least in 
the inner half of the radiative zone (compare the red and blue curves in Fig. H]), therefore 
such (enhanced) extra mixing has to be accompanied by an efficient heat transport. We 
have shown that the increase of F to a value between 0.75 and 1 also helps to explain the 
photometric peculiarities intrinsic to both the majority of red giants in tidally locked binaries, 
including the RS CVn stars, and the Li-rich K giants studied by KRL. However, post-bump 
luminosity stars with F ^ 1 would make zigzags that are too l e ngthy , in conflict with 
observations. This disagreement can be resolved if, following IZahru (119921 ). we assume that 
the enhanced extra mixing is produced by highly anisotropic turbulence whose associated 
chemical and heat transport in t he horizontal direction strongly domina tes over those in the 
vertical direction, D^ ^ -Dmix- IDenissenkov. fc Pinsonneaultl (l2008bf ) have demonstrated 
that, in this case, F = Dmi./[2a\K + Dh)] and / = 3T{Drai./K)/[l + T + 3T{D^i,/K)], 
while the expression V = (1 — /)Vrad+/Vad still holds. In the standard mixing length theory, 
F = 0.75 and F = 1 correspond to / = 0.66 and / = 0.75, respectively. For the anisotropic 
turbulent mixing, we obtain F ^ 1, assuming that Z^h ^ -Dmix ^ K and a > 1. At the 
same time, we can still get the observationally supported value of / ~ 0.75, provided that 
T{D,nijK) ^ 1, or 2{D^/K) ^ {D,^jKf. This means that for Drr.,^ ^ 1 0^ K (F ig. H, we 
need Dh ~ 5 x lO'^Dmix to keep / close to 0.75. Interestingly, iDenissenkovl ( 1201 Obi ) has used 
a similar ratio of D^ to -Dmix for the rotation-driven turbulent diffusion to make his model of 
magnetic braking of solar rotation consistent with observational data. Moreover, the same 
order of magnit ude ratios Dh I Dmn we re obtained for radiative zones of bump luminosity red 



giant models by iPalacios et al.l (|2006[ ) who took into account self-consistently the transport 
of angular momentum by rotation-driven meridional circulation and shear turbulence. 

Given that Vrad ~ 0.2 and Vad ~ 0.4 in the radiative zone of a bump-luminosity red 
giant, the maximum deviation of V from Vrad that can be achieved at the break-up rotation 
corresponds to / ~ 2ecrit = 0.48, which is too small to reproduce the extended zigzags 
presumably made by the Li-rich K giants and RS CVn giant components. However, this 
conclusion is based only on the analysis of the modification of the temperature gradient by 
the rotational distortion. It ignore s the fact that effects of rotation are also incorporated in 
other equations of stellar structure ( IDenissenkov. &: VandenBerg 112003b ). With all the effects 
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taken into account. iDenissenkov. Chaboyer. fc Li I ( 120061 ) were actually able to construct the 
evolutionary track of a 1.7 Mq star whose post-bump luminosity zigzag reached the HRD 
domain occupied by the RS CVn binaries. 

Comparing the green and red curves in Fig.HJ we infer that the ^He-driven thermohaline 
convection could be responsible for the observed Li enrichment only if the aspect ratio of 
its fluid elements were as large as a ~ 300. At present, we consider this highly improbable. 
The more likely interpretation from our point of view is that the Li enrichment in the 
stars studied by KRL and their location on the lower RGB below the bump luminosity are 
both caused by their fast internal rotation and its associated turbulent mixing and heat 
transport. This hypothesis is supported by the fact that surface rotation velocities of field 
stars with M > 1 . 6 Mc^ remain constant as they evolve on and away from the main sequence 
( jWolff. &: Simon 1119971 ). It is likely that only the most rapidly rotating stars, of which a 
handful have v sini > 200 kms~^, become Li-rich and make extended zigzags after the bump 
luminosity. Our hypothesis does not exclude the possibility of Li enrichment above the bump 
lumino sity in a red giant that has be en spun up as a result of its engulfing an orbiting giant 
planet (IDenissenkov. fc Weiss Il2000l ). 



The only possible scenario in which a substantial r aixing of H-rich material occur s dur- 
ing the He-core flash is the "hydrogen injection flash" (JMocak. Siess. &: Miiller Il201l[ ). The 
He ignition starts off-center, which leads to the formation of a He convective shell. For this 
scenario to work, convection from the He-burning shell must penetrate into the HBS. In the 
MESA stellar evolution code, such penetration, also known as convective overshooting, is 
modeled by the exponentially decaying diffusion coefficient Dqy = Dq exp[—2\r—ro\/{fHp)], 
where Hp is the pressure scale height, and Dq is a convective diffusion coeffici ent estimated 



using a mixing-length theory at the radius tq near the convective boundary (JHerwig et al. 



19971 ). Our computations of the He-core flash in the 1.7 Mq star show that H can be in- 
jected into the He convective shell only when / > 0.15. This exceeds the observationally 
constrained value of / by nearly an order of magnitude (JNoels et al.ll2010l . and references 
therein). Besides, given the discussed similarities between the Li-rich K giants studied by 
KRL and the RS CVn giant components and the fact that the latter stars are deflnitely 
on the lower RGB, because otherwise they would have undergone a common-envelope event 
with their close binary companions, we believe that these stars are in the same evolutionary 
phase. It should also be noted that the inclusion of convective overshooting in main-sequence 
stars would shift the maximum initial mass of the solar-composition stars that experience 
the He-core flash from M ^ 2.2 Mq to M ^ 1.8 Mq. 

Our hypothesis could eventually be tested by asteroseismology because high-precision 
photometry capable of measuring g-mode oscillations in red giants allows "to distinguish 



unambiguously between hydrogen-shell-burning stars (period spacing mostly ~ 50 seconds) 



and t hose that are also burning helium (period spacing ~ 100 to 300 seconds)" (JBedding et al. 



20111 ). 
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Fig. 1. — The black and green curves are the standard evolutionary tracks of the 1.7 Mq and 
2.2 Mq solar-composition stars. The red curves show the bifurcations caused by convective 
heat transport with the efficiencies F = 1 (the lower curve) and F = 0.75 (the upper curve) 
that starts in the radiative zone above the bump luminosity. As a result, the stars make 
excursions towards lower luminosities before resuming their ascents along the RGB. The 
blue and red star symbols represent, in turn, the Li-rich K giants studied by KRL and the 
giant components of RS CVn binaries. The neighbouring open symbols are separated by one 
million years. 
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Fig. 2. — The solid and dashed red curves show the changes of the surface Li abundance 
produced by enhanced extra mixing in the radiative zone of the 2.2 Mq star above the bump 



luminosity. It is assumed that D.^ 



2 X 10^1 cm2s-\ while V = 0.75 and T = 1 in 



the first and second cases. The green curve represents the evolutionary variations of the Li 
abundance caused by the FDU and thermohaline convection with the fluid element's aspect 
ratio a = 10. The blue star symbols are the Li- rich K giants. 
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Fig. 3. — The same as in Fig. |2l but for the surface carbon isotopic ratio. 
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Fig. 4. — The profiles of different diffusion coefficients in the radiative zone of the 2.2 Mq 
star at the bump luminosity. 



